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Abstract 

The relativistic "no pair" model of quantum electrodynamics uses the Dirac operator, D{A) 
for the electron dynamics together with the usual self-energy of the quantized ultraviolet cutoff 
electromagnetic field A — in the Coulomb gauge. There are no positrons because the electron 
wave functions are constrained to lie in the positive spectral subspace of some Dirac operator, 
D, but the model is defined for any number, N, of electrons, and hence describes a true many- 
body system. In addition to the electrons there are a number, K , of fixed nuclei with charges 
< Z. If the fields are not quantized but are classical, it was shown earlier that such a model is 
always unstable (the ground state energy E — -co) if one uses the customary D(0) to define 
the electron space, but is stable (E > — const. (N + K)) if one uses D(A) itself (provided the fine 
structure constant a and Z are not too large). This result is extended to quantized fields here, 
and stability is proved for a = 1/137 and Z < 42. This formulation of QED is somewhat unusual 
because it means that the electron Hilbert space is inextricably linked to the photon Fock space. 
But such a linkage appears to better describe the real world of photons and electrons. 

1 Introduction 

The theory of the ground state of matter interacting with Coulomb forces and with the magnetic 
field is not yet in a completely satisfactory state. Open problems remain, such as the inclusion of 
relativistic mechanics into the many-body formalism and the inclusion of the self-energy effects of 
the radiation field, especially the quantized radiation field. 
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One of the fundamental attributes of quantum mechanics is the existence of a Hamiltonian with 
a lowest, or ground state energy, and not merely the existence of a critical point of a Lagrangian. 
The 'stability' problem, which concerns us here, is to show that the ground state energy is bounded 
below by a constant times the total number of particles, N + K, where N is the number of electrons 
and K is the number of nuclei - whose locations, in this model, are fixed, but chosen to minimize 
the energy. We do not discuss the existence of a normalizable ground state eigenfunction, as in || , 
but only the lower boundedness of the Hamiltonian. 

This problem has been resolved successfully in various models such as the usual nonrelativistic 
Schrodinger Hamiltonian with only electrostatic interactions. Further developments include ex- 
tensions to relativistic kinetic energy \/p 2 + m 2 — m in place of the nonrelativistic p 2 /2m, and 
extensions to matter interacting with classical magnetic fields (including a spin-field interaction 
B), stabilized by the classical field energy 



^ J B(x) 2 dx , (1) 



and then the quantization of the B field. Many people participated in this development and we 
refer the reader to ||T7j| and the references therein for an account up to 1997. 

In this paper we take a further step by addressing the problem of relativistic matter, using the 
Dime operator (without pair production, i.e., the "no-pair" model) interacting with the quantized 
radiation field having an ultraviolet cutoff A. In (I?]] the corresponding problem was solved with 
a classical radiation field, in which the field energy is given by (H), and we shall use some of the 



ideas of that paper here. The idea for such a model goes back to and [23]. With a classical B 
field no ultraviolet cutoff is needed, but it is needed with a quantized field, for otherwise the field 
energy diverges. 

Because of the ultraviolet cutoff our model, which in other respects is relativistic, is not truly 
relativistic at energies of the order of the cutoff. We have not, however, attempted to renormalize 
the theory and, if this can be done consistently, the resulting theory will be relativistic at all 
energies. 

In Q the problem of nonrelativistic electrons (with spin) interacting with the quantized ultra- 
violet cutoff field was solved by using results in |l(| but using only the part of the field energy 
within a distance 1/A of the fixed nuclei. The constants and exponents in Q were improved in 
in particular, the Hamiltonian is bounded below by — AK. The relation of the classical field energy 
to the quantized field energy involves a commutator that, when integrated over the whole space R 3 
yields an infinite constant, even with an ultraviolet cutoff. This is the reason for considering only 
a local field energy, since only a local field energy yields a finite commutator, and we do the same 
here. 



In Section |2| our model is defined and the main Theorem 2A is stated. With the fine structure 



constant a = 1/137, stability holds for Z < 42. The main idea of the "no-pair" model is that 
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there are no positrons, and electronic wave functions are allowed to lie only in the positive spectral 
subspace of some Dirac operator D. While the Dirac operator D(A), which is contained in the 
Hamiltonian and which defines the electron dynamics, always contains the magnetic vector potential 
A(x), the operator D that defines an electron could be -D(O), the free Dirac operator. Indeed, this 
is the conventional choice, but it is not gauge invariant and always leads to instability as first shown 
in [17] for classical fields and here for quantized fields. 

The question of instability is complicated. There are two kinds (first and second) and two cases 
to consider (with and without Coulomb potentials). Instability of the first kind means that the 
ground state energy (bottom of the spectrum of the Hamiltonian) is — oo. Instability of the second 
kind means that the energy is finite but is not bounded below by a constant times N + K. The 
occurence of these instabilities may or may not depend on a and Z and whether or not a cutoff A 
is present. 

The physical nature of the instability, if it occurs, is different in the two cases. When it occurs 
in the absence of Coulomb potentials (meaning that the aV c term in (|ll]) is omitted) it is due to 
the y/aA(x) term in D{A) blowing up. When it occurs because of the Coulomb potentials being 
present it is due to an electron falling into the Coulomb singularity of the nucleus. The various 
possibilities, all proved in this paper, are summarized in detail in the following two tables and 
discussed in detail in Appendix ||. For the proofs of the instabilities listed here, we rely heavily on 
and @. 



Electrons defined by projection onto the positive 
subspace of -D(O), the free Dirac operator 





Classical or quantized field 


Classical or quantized field 




without cutoff A 


with cutoff A 




a > but arbitrarily small. 


a > but arbitrarily small. 


Without Coulomb 


Instability of 


Instability of 


potential aV c 


the first kind 


the second kind 


With Coulomb 


Instability of 


Instability of 


potential aV c 


the first kind 


the second kind 
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Electrons denned by projection onto the positive 
subspace of D(A), the Dirac operator with field 





Classical field with or without cutoff A 
or quantized field with cutoff A 


Without Coulomb 
potential aV c 


The Hamiltonian is positive 


With Coulomb 


Instability of the first kind when either 
a or Za is too large 


potential aV c 


Stability of the second kind when 
both a and Za are small enough 



The main point of this paper is the proof of the bottom row of the second table in the quantized 
case; the classical case was done in [17]. 

There are several ways in which one could hope to go further. One is that one should really 
prove stability for the binding energy, i.e., one should compute the energy difference between that 
of free particles and that of the interacting system. In a theory with quantized fields the self-energy, 
i.e., the energy of a free electron, is unknown and quite large. As we show in jOJ and Jlij the self- 
energy of a nonrelativistic particle with spin is bounded below by +A, and probably even +A 3 / 2 . 
Moreover, for N fermions (but not for bosons) this energy is proportional to C'NA with C > 0. 

Another very important problem to consider is renormalization; our mass m is the unrenormal- 
ized one. An answer to this problem also has to address the question of the meaning of mass in 
an ultraviolet cut-off model, since several definitions are possible. Is it the coefficient of j3 in an 
effective Dirac operator that gives the renormalized dynamics, or is it the ground state energy of a 
"free" electron? 



The results in this paper will be used in another paper of ours |L5| to give upper bounds to the 
hydrogen atom binding energy (and hence to the mass renormalization using the first definition) 
in this relativisitic no-pair model and in some non-relativistic models with quantized fields. 

Finally, let us note that the inclusion of positrons into the model cannot change the fact that 
defining an electron by means of D(0) will still cause the instabilities listed in the tables above. 
The reason is simply that the existence of positrons does not prevent one from considering states 
consisting purely of electrons, and these alone can cause the listed instabilities. 

The use of D(A) instead of D(0) to define the electron requires a significant change in the 
Hilbert space structure of QED. It is no longer possible to separate the Hilbert space for the 
electron coordinates from the Hilbert space (Fock space) of the photons. The two are now linked 
in a manner that we describe in the next section. 
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2 Basic Definitions 

We consider N relativistic electrons in the field of K nuclei, fixed at the positions R±, Rk £ K 3 - 
(In the real world the fixed nuclei approximation is a good one since the masses of the nuclei are 
so large compared to the electron's mass.) We assume that their atomic numbers Z\, Zk are all 
less than some fixed number Z > 0. Since the energy is a concave function of each Zj separately, 
it suffices, for finding a lower bound, either to put Zj = 0, i.e., to remove the j-th nucleus, or to 
put Zj = Z (see ||). Thus, without loss of generality, we may assume that all the nuclear charges 
are equal to Z. 

We use units in which h = 1 and c = 1. a = e 2 /he is the dimensionless "fine structure constant" 
(=1/137 in nature). The electric charge of the electron in these units is e = \fa. 

We use the Coulomb, or radiation gauge so that the Coulomb potential is a function only of 
the coordinates of the N electrons, xi, X2, ... ,xjsj and equals aV c , where 

^EE^-^ E E w 

i=l k=l 1 K ' l<i<j<N 1 Jl l<k<l<K 1 K 11 

In this gauge, it is the vector potential that is quantized. A careful discussion of the field and its 
quantization is given in Appendix [A]. The (ultraviolet cutoff) magnetic vector potential is defined 
by 

A(x) = -L J2 [ 4= U(k)e ik -* + aj(fc)e-**) dk , (3) 

27r J\k\<A Joj(k) v / 



where A is the ultraviolet cutoff on the wave-numbers \k\. The operators a\,a* x satisfy the usual 
commutation relations 

[a x {k),a* u (q)} = 5(k - q)5 x ,u , [a\(k), a u (q)] = 0, etc (4) 

and the vectors e\(k) are the two possible orthonormal polarization vectors perpendicular to k and 
to each other. 

Our results hold for all finite A. The details of the cutoff in @ are quite unimportant, except 
for the requirement that rotation symmetry in fc-space is maintained. E.g., a Gaussian cutoff can 
be used instead of our sharp cutoff. We avoid unnecessary generalizations. The cutoff resides in 
the A-field, not in the field energy, Hf, sometimes called c£T(o>), which is given by 

H f = E / u(k)a* x (k)a x (k)dk (5) 

A=l,2 

The energy of a photon is uo(k) and the physical value of interest to us, which will be used in 
the rest of this paper, is 

u(k) = \k\ (6) 
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Again, generalizations are possible, but we omit them. 

An important fact for our construction of the physical Hilbert space of our model is that 
[A(x),A(y)] = [B(x),B(y)] = [A(x),B(y)] = for all x,y. Here, B is the magnetic field given by 

2 

B(x) = curl A(x) = j-i2f ~TlTr i a ^y k ' X ~ *l(k)e- lk - x ) dk . (7) 

The kinetic energy of an electron is defined in terms of a Dirac operator with the vector potential 
A{x) (with x being the electron's coordinate) 

D(A) := ct ■ (-iV + ^A(x)) + m/3 , (8) 

with a and (3 given by the 2x2 Pauli matrices and 2x2 identity I as 



Note that 

D{A) 2 = f p {A){A)+m 2 , (9) 



where T P (A) = I ^ 1 tP J and T P (A) is the Pauli operator on L 2 (R 3 ;C 



2\ 



T P {A) = [a- (p + ^A(x))] 2 = (p + ^tA(x)f + ^a-B(x) . (10) 

As a step towards defining a physical Hamiltonian for our system of N electrons and K fixed 
nuclei, we first define a conventional, but fictitious Hamiltonian 

N 

H' N = Y,Di(A)+aV c + H f . (11) 
i=i 

This H' N acts on the usual Hilbert space Hn = &) N L 2 (R 3 ; C 4 ) ® J 7 , where T is the Fock space 
for the ^4-field. A vector in Hn can be written 

CO 

* = (J)$j(xi, ...,x N ; T!,...,t n ; ki,...,kj] Xi,...,Xj) . (12) 

3=0 

Here, the Aj take the two values 1,2 and the the n take the four values 1,2,3,4. Each is 
symmetric in the pairs of variables fcj,Aj and it is square integrable in x,k. The sum of these 
integrals (summed over A's, r's, and j) is finite. The operators a\(k) and their adjoints act, as 
usual, by 

oo 

a\(k)^ = (J) y/j + l$ J+ i(xi, ...,xat; ri,...,rjv; h, ...,kj,k; X±,...,Xj,X) . (13) 

j=0 
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As mentioned before, the physical Hilbert space is constructed using the positive spectral pro- 



jections of the Dirac operators D{(A). By Lemma C,l the N Dirac operators commute in the strong 
sense that their spectral projections commute with each other. Thus, the Hilbert space TCn can be 
divided into 2^ subspaces according to the positive and negative spectral subspaces of each Di(A). 
(Note that as long as m > there is no zero spectral subspace.) We denote by P + the orthogonal 
projection onto the positive spectral subspace for all the Dirac operators. 

The space P + T~Ln is invariant (up to unitary equivalence) by the natural action of the permu- 
tation group Sjsr consisting of permutations of the electron labels. In accordance with the Pauli 
principle we choose the antisymmetric component of P + TCn, as the physical Hilbert space. Thus, 
our physical Hilbert space is given as 

n fys =AP + nN (14) 

where A is the projector onto the antisymmetric component. 

Formally, i.e., without attention to domain questions, our physical Hamiltonian on 7i^ ys is 
defined to be 

^phys = p+ H > N p+ . ( 15 ) 

Since we are interested in this operator as a quadratic form, it suffices to specify a domain 
Qn which is dense in 7i v ^ ys and on which the expectation values of the all the operators involved 
are finite. Since all the operators are symmetric, and since a stability estimate entails that the 
quadratic form is bounded below, its closure exists and defines a selfadjoint operator H^ ys . Such 
a domain Qn is constructed in Appendix |D|. Note that by definition Q consists of antisymmetric 
elements. 

We note that each of the ./V Dirac operators commute with P + . For ip £ Qn we have Di(A)ip = 
P + Di(A)ip. For the other two terms in ( |i~l| ) the role of the projector is not so trivial and that is 
why we have to write P + H' N P + . 

This model has its origins in the work of Brown and Ravenhall pi] and Sucher p3|]. The 



immediate antecedent is [17]. 
Let us note five things: 

(i) . It is not entirely easy to think about 7i^ ys because the electronic L 2 -spaces and the Fock 
space are now linked together. In our choice of positive energy states, the electrons have their own 
photon cloud. We chose to apply the projector P + first and then antisymmetrize. As explained 
in Appendix O we can, of course, do it the other way around and obtain the same Hilbert space, 
since P + commutes with permutations. We also show in Appendix ^ that 7i^ ys is not trivial; in 
fact it is infinite dimensional. 

(ii) . Usually, in quantum electrodynamics, one defines TL^ ys by means of the positive spectral 
subspace of the free Dirac operator D(0) = —ia ■ V + mf3, instead of D(A). This is easier to think 
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about but, as demonstrated in [|17]] with a classical A field instead of a quantized field, the choice 
of D(Q) always leads to instability, as listed in the tables in Section || and discussed in detail in 
Appendix |E|. 

(hi). Because of the restriction to the positive spectral subspace of D(A), the Dirac operator 
is never negative. The only negative terms in H^ ys come from the Coulomb potential. It should 
also be noted that the choice of the free Dirac operator to define an electronic wave function is not 
a gauge covariant notion. The D(A) choice is gauge covariant. 

(iv) 7i^ ys depends on a and m. 

(v) While energy, being one component of a four-vector, is not a relativistically invariant quan- 
tity, it is true, nevertheless, that positive and negative energies of D(A) are relativistic concepts 
since they are invariant under Lorentz transformations that do not change the direction of time. 
We thank J-M. Graf for this remark and we thank J. Yngvason for noting that for this to be true 
it is essential that the joint spectrum of energy and momentum of D(A) lies in the light cone. We 
have not proved this, but it is plausibly true. 

is 



Our main result, to be proved in Section 2.1 



2.1. THEOREM (Relativistic Quantum electro dynamic Stability). Assume that Z and 



a are such that there is a solution k and e > to the three inequalities (pq), ( pq j and (5^). 
Then Hfj iys in ftT^) is bounded below by 

H p N hys >+V^mN- ^KCl , (16) 

where 

^_N Qy/T^~e + (q/2)(vgg + 2.3) 2 

° 2 " K 27^ • U?j 

In particular, Z < 42 is allowed when a = 1/137. 



Actually, our proof of Theorem 2.1 utilizes the absolute value of the Dirac operator |D(A)| on 
the Hilbert space AH.N- If we recall the connection between D{A) and the Pauli operator in (^) 
and (0), we can prove the following theorem as a byproduct of our proof of Theorem |2.1| . Note 
C 2 in place of C 4 here. 

2.2. THEOREM (Stability with the Pauli operator). Let H'^ = V tP ( A ) + ™ 2 +<xV c + 

Hf be a Hamiltonian on the space A &) N L 2 (M 3 ; C 2 ) (^) T . Then stability of the second kind holds 
under the conditions stated in Theorem [g. jL and with the same lower bound ([7^). 



3 Bounding the Coulomb Potential by a Localized Relativistic 
Kinetic Energy 



The following Theorem 3^ contains the main technical estimate needed in this paper, but it is 
independently interesting. It deals with a model of relativistic electrons interacting with quantized 
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fields, but without the spin-field interaction and without the field energy. While this model is 
different from the no-pair Hamiltonian (|l5|), some of its properties will be useful later. We consider 
two such Hamiltonians: A usual one 

N 

K^2\Pi+ VaA(xi)\ +V C , (18) 
i=l 

(with k > 0) and a related one with a localized kinetic energy described below in (|24|), (p5|). In this 
section A{x) is some given classical field, not necessarily divergence free. There is no a in front of 
V in (|l8|). The Hilbert space is A (gj^L^IR 3 ;^) for fermions with q 1 spin states'. 

With the K nuclei positioned at distinct points Rj G M 3 , for j = 1,...,K, we define the 
corresponding Voronoi cells by 

Fj = {xeR 3 : \x-Rj\ < \x-Ri\,i = l,...,K,i^j} . (19) 

These Voronoi cells are open convex sets. We choose some L > and define the balls Bj C M 3 by 

Bj = {x : \x-Rj\ < 3L} , (20) 

and denote by B the union of these K balls and by Xb the characteristic function of B. Similarly, we 
define smaller balls, Sj = {x : \x — Rj\ < 2L}, and define Xs to be the characteristic function of the 
union of these K smaller balls. Choose some function g £ VF 1 ' 1 (R 3 ) with support in {x : \x\ < 1}, 
with g > and with f g = 1. Define gz(x) = L~ 3 g(x/L). Clearly J gi, = 1 and gL has support in 
{x : |x| < L}. With * denoting convolution, set 

Mx) =9L*Xs{x) ■ (21) 

This function (pi is nonnegative and everywhere bounded by 1. We also define 4>2 = 1 — 4>i and set 

F = <t>i/yJ + and G = H/^l + 4>\ ■ (22) 

Note that (f>\ + </>| > 1/2 and F(x) = 1 if |x — Rj\ < L for some j. Note also that F and 4>i are 
supported in B, i.e., Xe^i = 0i and XbF = F- 
We find that 

|VF| 2 + |VG| 2 <4|V^| 2 <-^^ a |V 5 (x)|^ . (23) 

and hence |VF|,|VG| < 2| V0i|. 

The function g that minimizes the integral in (p3|) is g(x) = 3/4ir for \x\ < 1 and zero otherwise. 
(Although this <? is not in VK 1 ' 1 (M 3 ) it is a limit of T^ 1,1 (M 3 ) functions.) Then the integral equals 
3 and |VF| 2 + |VG| 2 < 36L~ 2 . 

The localized kinetic energy operator Q(A) is given by 

Q(A) = F(x) \p + y/aA(x)\ F(x) = F(x)\J {p + yfaA(x)) 2 F{x) . (24) 
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This operator is well defined as a quadratic form since the function F is smooth, and hence defines 
a self adjoint operator via the Friedrichs extension. 

The related relativistic Hamiltonian, with localized kinetic energy, is now defined by 

N 

H% c :=kJ2QM) + V c , (25) 
i=l 

and has the following bound which, it is to be noted, does not depend on the details of g(x). 

3.1. THEOREM (Bound on Coulomb energy). For any vector field A{x) and for N fermions 
with q spin states, 

N 

OJA) + V„ > - .,„ 

2L LV ' ' 21 1 ~ 2L 



K^Q J (A) + y c >-^max{(V2Z + l) 2 , 2Z + ^} > -^(v^Z + 2.3) 2 , (26) 

8=1 

provided k > max{g/0.031, irZ}. 

Proof. It was proved in (eqns. (2.4-2.6) with A = 10/11) that the Coulomb potential V c is 
bounded below by a single-particle potential plus a constant, namely, for Xj, Rj S M 3 , 

TV ^ 2 K 



iV Z 2 K 



i=i ' j=i 3 



where 2Dj = minj^{|xj — Xj\} and, for x £ Tj, 



= (^+iA/2) 2 !oilx _ Rl> m 

y ' \x-Rj\ 3 ~ 11 

Z 121 „ , . WD, , . 

+ 77777" for\x-R j \<— ±. (28) 



\x-Rj\ A2Dj 1 -"1.1 



This estimate reduces our problem to finding a lower bound to 

N N N y 2 K i 

K £ Q,(A) - £ F( Xi ) 2 W( Xi ) - £(1 - FCxi) 2 )^) + • 

i=l i=l i=l j=l 3 

Since F(x) = 1 if \x — Rj\ < L for some j, the third term in (^) is bounded below by 

N r , nrrz . ^ 2 ~~ 110 - 



2L max{(V2Z + l) 2 , 2Z + — } . (30) 

Estimating the first and second terms using the Pauli exclusion principle amounts to filling the 
lowest possible energy levels with q electrons each, and this energy is bounded below by q times 
the sum of the negative eigenvalues of the operator 

F{x)(\p + y/aA(x)\ - W{x))F{x) . (31) 



According to the generalized min-max principle |12| Corollary 12.2, and the fact that H-Fi/'ll < 
, this is bounded below by q times the sum of the negative eigenvalues of the operator | (p + 
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\faA{x))\— W{x). However, Theorem 1 of [19], shows that this sum is not less than (-Z 2 /8) £* =1 1/Dj 
under the stated condition on k. 

(Notes: We refer here to Theorem 1 of [19] because, as noted in [17], the proof of that theorem 
holds for \p + ^/aA{x)\ in place of \p\. While Theorem 1 of |jl]] is stated in terms of V c , the proof 
in [19] actually replaces V c by its lower bound pTj), ) □ 



4 Proof of Theorem [27T| 

We employ a strategy similar to that in 



As a first step we use Theorem 3.1 with a suitable choice of L to control the Coulomb potential. 



The operators appearing in Theorem 3J do not involve spin, but the number of spin states, 
q, is important for determining the relevant value of k. The correct choice is q = 2, not q = 4, 
as explained in |17j page 42 and appendix B. The point is the following. The one-body density 
matrix T(x, a; x', a') coming from an antisymmetric N particle wave function ^ defines a reduced 
one body density matrix 

4 

7 (x,x') = ^r(x,cr;x» . (32) 

<T=1 

This reduced density matrix, in general, satisfies < Tr7 = J "f(x,x)dx < 4. If, however, ^ is 
in the range of P + , then < Tr7 < 2, as shown in [17]. In the proof of Theorem 3T, the only 
relevant information about ^ enters via the reduced single particle matrix 7. Thus, we require only 
K > max{64.5, ttZ}. 

In the definition of F we set L = C2/A where C2 > is some constant to be conveniently chosen 
later. We then have (recalling (§), (M), and P l + D i (A)P l + = P+| A(-4)|^ + ) 



A' 



Y^Di(A) + aV c 



,i=l 



iV 



i=i 



Tf{A)+m 2 -KaQi{A) 



aN^iV^Z + 2.3) 2 P+ . 
2G2 

(33) 



(Here, Q(A) really denotes the 4x4 operator Q(A) <g> I4 where J4 is the identity in spin-space.) 
Consider the operator 

N 



H 2 :=P+Y, 



TP {A) +m 2 -5m- KaQi(A) + C 3 A 



P + 



(34) 



where the numbers < 5 < 1 and C3 > will be chosen later. 

If we denote by H + the projection onto the positive spectral subspace of D{A) acting on 
L 2 (M 3 ;C 4 ) (g)^", then H 2 is bounded below by 



Tr 4 E+,SH H 



(35) 
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where Tr n with n = 1, 2,4 denotes the trace on L (R ; C n ). The operator S is 



S ■= JT P (A) + m 2 -5m- KaQi(A) + C 3 A . 



It has the form 



S 



'Y N 
.0 Y . 



(36) 



(37) 



Here, the entry Y is a 2 x 2 matrix valued operator and [AT]_ denotes the negative part of a self- 
adjoint operator X (and which is nonnegative by definition). The projection E + is not explicitly 
given, but observing, as in [T^j , that the projection H~ onto the negative energy states is related to 
~+ by 



U-t-E+U = -UE + U , 



where U is the matrix 



U 



o r 

-I o, 



we see that the operators and E SE have the same spectrum. Thus, 

1. 



Tr 4 [H+SH+]_ < -Tr 4 [5]_ =Tr 2 



T P (A) +m 2 -5m- kcxQ{A) + C 3 A 



(38) 



(39) 



(40) 



Therefore, the infimum of the spectrum of H% over states that satisfy the Pauli exclusion principle 
(with 4 spin states) is bounded below by 



-Tr 2 



T P (A) +m 2 -5m- KaQ(A) + C 3 A 



(41) 



The BKS inequality Q (see also [[[7]]) states that for positive operators A and B, Tr 2 [A — £?]_ < 
Tr 2 [vl 2 - S 2 ]i /2 . Note that y/T p (A) + m 2 - 5m > and, therefore, 

1/2 



ffo > -Tr 2 



T P (A) + m 2 - dm + C 3 A - k a Q(A) 



(42) 



which is greater than 



-Tr 2 



T p (A) + m 2 - 5m^) + C 3 2 A 2 - K 2 c?Q{Af 



1/2 



(43) 



1/2 

(Here, and in the following, we use the fact that Tr[X]_ is monotone decreasing in X.) 

Next, we expand (• • • ) 2 in ( f43| ) and use the arithmetic-geometric mean inequality to bound 
from below by 

-Tr 2 [(T P (A) + m 2 )(l - e) + (1 - l/e)m 2 5 2 + C 2 A 2 - K 2 a 2 Q(A) 2 ] 1/2 . 



(44) 
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We choose 5 so that the mass disappears, i.e., 5 2 = e. 

The next step is to localize the Pauli term T P (A). A standard calculation shows that (with 
F, G as in Section ||) 

T P {A) = FT P (A)F + GT P {A)G - |VF| 2 - |VG| 2 > FT P {A)F - \VF\ 2 - \VG\ 2 . (45) 
We insert the right side of (|45|) into ( EjJ ) and, recalling (p3|), choose C3 to eliminate the A 2 term, 



i.e., 



Cs = ^ / |V 9 (x)|dx . (46) 

Thus, using the fact that Q(^) 2 = ^/aA(x)|F 2 |p + ^/a^(x)|F < F(p+ *faA(x)) 2 F, which 

follows from F 2 < 1, we obtain the bound 

H 2 > -Tr 2 [(1 - e)FT p (A)F - K 2 a 2 Q(A) 2 ] V ' 

> -Tr 2 lF{{l-e-K 2 u 2 )(p + ^A{x)) 2 + {l-e)^XBO-B{x))F] V2 . (47) 

We have used the fact that X13F = F. 

Since FXF > — F[X]-F for any X, the eigenvalues of FXF are bounded below by the eigen- 

1/2 1 /2 

values of — F[X]_F, and hence we have that Tr [FXF]_ < Tr [FX^F] ' , and hence 
Tr 2 [F{(l-e- K 2 a 2 )(p + ^A{x)) 2 + (1 - e)JaXB° • ^] V2 

1 /2 

< Tr 2 {F [(1 - e - K 2 a 2 )(p + ^A(x)) 2 + (1 - e)^X B a • B(xj\ _ f} . 

The expression [ ]_ between the two F's is, by definition, a positive-semidefinite self-adjoint 
operator and we denote it by Y. Now 

Tt 2 {FYF) 1 / 2 = Ti 2 {FY l / 2 Y 1 / 2 Ff/ 2 = Ti 2 (Y 1 / 2 FFY 1 / 2 ) 1 / 2 (48) 

since, quite generally, X*X and XX* have the same spectrum (up to zero eigenvalues, which are not 
counted here). Finally, we note that since F 2 < 1, Y l / 2 FFY 1 / 2 < Y, and hence Ti 2 (FYF) 1 / 2 = 
Tr 2 (y 1 / 2 i ? i ? y 1 / 2 ) 1 / 2 < Tr 2 y x / 2 . Thus, it remains to find an upper bound to [h]]/ 2 where 

h = (1 — e — K 2 a 2 ){p+^A{x)) 2 + {1-e)^XbvB{x) . (49) 

Denote the negative eigenvalues of h by —e\ < — e 2 < • • • . One way to bound the eigenvalues 
from below is to replace a ■ B(x) by — \B(x)\, but then each eigenvalue of h := (1 — e — n 2 a 2 )(p + 
T/aA(x)) 2 — (1 — e)y/aXB\B{x)\ on L 2 (M 3 ) would have to be counted twice (because Tr 2 is over 
L 2 (R 3 ;C 2 ) and not L 2 (R 3 )). As shown in j^0|, however, the intuition that each negative eigenvalue 

' — ' 1/2 — 1/2 

of h should be counted only once is correct. Thus, Tr 2 [/i]_ < Tri[/t]_ . By the Lieb-Thirring 
inequality we obtain the bound 

E < n (1 " g) ?" /2 / . ( 5 °) 

~ (1 — e — k or ) 1 Jb 
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with £ = 0.060 |0|. 

It is to be emphasized that (50) is an operator inequality. That is, the operator in (|34|), which 
is part of Hf^ ys , satisfies 



Ho > 



(1 -e) 2 I a 
(1 - e - K 2 a 2 ) 3 / 2 



B(x) 2 dx . 



(51) 



The right side of ([5l]) can be controlled by the field energy through inequality ( |79| ) — provided 
1/87T is not less than the constant in (|50|), (|5l|). 



4.1 Evaluation of Constants 

We are now ready to list the conditions on the constants C2 and e that have been introduced and 



to use these to verify the results of Theorem 2.1. 



k > max{64.5, trZ} (52) 

Conditions : (na) 2 < 1 — e < 1 (53) 

(l-e) 2 q 1 , 41 

(1 -£- k 2 o 2 ) 3 / 2 ~ 8^(0.060) 1 ' 



The first comes from Theorem 3.1 with q = 2. The second is the condition that the kinetic energy 
term in H2 is positive. The third is the requirement that the the field energy Hf dominates the 
sum of the negative eigenvalues in (|50|) . 

Assuming these conditions are satisfied the total energy is then bounded below by the sum of 
the following four terms (recalling (03) and 5 2 = e): 



emN (55) 
^AiV (50) 



Energy Lower Bounds : — — ( v2Z + 2.3 ) N (57) 

2C/2 



~ I 1 > -- A (3C 2 ) 3 K = --AC! K (58) 
8vr 2 J B ~ 3 8vr 2 v ; 2vr 2 v ; 

The first comes from the —5m term in (|34"|). Similarly, the second comes from the +C3A term in 



(|34|). The third term is the last term in ( J33| ) which, in turn, comes from Theorem 3.1. The fourth 
term is the additive constant in (|79|) with w(y) = Xb{v)- The volume of B is bounded by the 
number of nuclei times the volume of one ball of radius 3L around each nucleus. 
Obviously we choose 

nA _ N 6VT^ + (a/2)(y / 2Z + 2.3) 2 

° 2 ~ K 27/2^ • (59) 
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The sum of the terms (55 - 58) then become our lower bound for the energy 



E _ 18A K „ 

N>+^m-—-Ci, (60) 

which satisfies stability of the second kind. 

To find the largest possible Z for which stability holds we take a = 1/137 and make the choice 
e = 0. We then find, from (|54j), that na < 0.97. Setting k = irZ we find stability up to Z = 42. 

The choice e = makes the energy in (^0[) negative. Recall that if Z = then E/N = m. To 
make contact with physics we would like the energy to be positive, i.e., only a little less than Nm. 
To fix ideas, let us consider the case irZ < 64.5 and a = 1/137. Then k = 64.5, k« = 0.471 and 
(kq) 2 = 0.222. From (0), we require that (with x = l- e> 0.222) 

x 2 < 90.9(x - 0.222) 3/2 , (61) 

which means that we can take 1 — s = .229 or e = 0.771. 

Now let us consider the case of hydrogen, Z = 1 and N = K (neutrality). From (^) we find 
that C 2 = 0.908. Then @ becomes 

E 

— = 0.866m - 4.29A . (62) 

If A is less than one fifth of the electrons's self-energy, the total energy of arbitrarily many hydrogen 
atoms is positive. This bound could be significantly improved by more careful attention to our 
various inequalities. 

A Appendix: A Note About Units 

The choice of units in electrodynamics is always confusing, especially when interactions with charged 
particles are involved. 

The interaction of the magnetic vector potential with a charged particle is eA(x). In cgs units 
the classical field energy is 



^classical = J_ f { B {xf + E(xf\ dx 

87r y R3 



(63) 



With B{x) = curL4(x), we use the Coulomb (or radiation gauge) so that div^4(x) = and divE'(x) = 
0. 

We define a\{k) and its complex-conjugate (classically) or adjoint (quantum-mechanically), 
a* x (k), in terms of the Fourier transform of (the real fields) A(x) and E(x) as follows. 

M*) = ^ (a X (k)e ik x + al(k)e-^) dk (64) 

E{x) = VW\ex(k)(a x (k)e ik - x -al(k)e- lk - x )dk, (65) 
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in terms of which 

B(x) = ^pE / ^ff 1 U(k)e ik -* - aJ(fc)e-*») rffc . (66) 

The parameter ^/Hc/2tt in (|64|-|66|) were chosen purely for convenience later on. The two unit vectors 
here, £\(k), A = 1,2, are perpendicular to each other and to k (which guarantees that div^4 = 0). 
They cannot be defined on the whole of M 3 as smooth functions of k (although they can be so 
defined with the use of 'charts'), but that will be of no concern to us. 

Thus, when - |66| ) are substituted in ( |6"3| ) we obtain (using Parseval's theorem and 
/ exp(ik ■ x)dx = (2ir) 3 5(k) and \k\ 2 e x (k) ■ s x (-k) = - (k A e x (k)) ■ (-k A E\(—k))) 

1 2 r 

^classical = V / \k\ {a* x (k)a x (k) + a x {k)a\ (k)} dk (67) 

(Although a x (k)a x (k) = a x (k)a* x (k) for functions, this will not be so when a x (k) is an operator. 
The form in (^) is that obtained after the substitution just mentioned.) 

To complete the picture, we quantize the fields by making the a x (k) into operators with the 
following commutation relations. 

[a x (k), a*y(k')] = 5 x , x ,5(k - k') and [a x (k), a y (k')] = . (68) 

The quantized field energy is obtained from (^37|, |6^) and is given by the Hamiltonian operator 

2 r 

H f = h c y" \k\a* x (k)a x (k)dk . (69) 

It agrees with ( |67|) up to an additive 'infinite constant'. 

In the rest of this paper we omit he since we use units in which h = c = 1. 



B Appendix: Field Energy Bound 



In this appendix we prove (|7y) which relates the localized classical field energy to the quantized 
field energy. A proof was given in Q . The small generalization given here is a slightly modified 
version of that in |l3|, 14 1 . 



Consider a collection of operators (field modes), parametrized by y G M 3 , and by j in some set 
of integers (j G {1,2,3} in our case of interest) given, formally, by 



2 f 

L M = E / VW\vxj(k)e ik -y ax (k)dk 

x=1 J\k\<A 



(70) 
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where v\ j is the Fourier transform of some arbitrary complex function v\ j(x). Our convention for 
the Fourier transform of a general function g(x) is 

g(k) = (2tt)- 3/2 [ g{x)e +ik - x dx and g{x) = (2tt)- 3/2 [ g{k)e~ lk - x dk . (71) 



The following lemma is elementary. It involves v X j(x) and a summable function w(x), with a 
norm defined by 

„ „ Ej fRz\Exf\*v\,j(x)\ 2 \w{x)\dx 

\H v ■= sup , (72) 

h(x) l^xJ^\fx(x)\ 2 dx 

where * is convolution. 



B.l. LEMMA (Lower bound on field energy). Assume that \\w\\ v < 1. Then 

H f >T [ w{y)L*(y)L 3 {y)dy . 



(73) 



j 

Moreover, if w{y) > 0, for all y then 



2 

H f^lT, ± [ w(y)(L j (y)±L*(y)) 2 dy-l^2fl I \ k \ \^M?dk [ w{y)dy , (74) 

4 , JR 3 Z ,T~ l J\k\<A JM 3 



3 A=l 



/or any choice of + or - for each j. (Note that —(L — L*) >0.) 



Proof. The difference of the two sides in (|73D is a quadratic form of the type 

X^a A' / / a x(k)Q(k, A : k' , \')a x ,(k')dkdk' . In order to establish ([73]) it is necessary and sufficient 
to prove that the matrix Q(k, A : k', A') is positive semidefinite. This is the condition that 

V/ \f x (k)\ 2 dk>(2Kf/ 2 y] E // fx(k)vxAk)fy(k')vyAk')w(k' - k)dkdk' 

(75) 

for all L 2 functions fx(k). Condition ( |75| ) is just the condition that \\w\\ v < 1, since / * v(k) = 
(2nf/ 2 f(k)v(k). 

To obtain from (73) we use the three facts that u;(x) > 0, that 

L*(y)L,(y) = Lj(y)Lj(y) - £ / |fc||«A,#)| 2 dA; , (76) 

J\k\<A 



A=l 

and that, quite generally for operators, 



± (LL + L*L*) < L*L + LL* . (77) 

□ 
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The following examples are important. First, we define the ultraviolet cutoff fields A A , B A , E A 



as in (p4|,p€,65) except that the k integration is over |fe| < A instead of R . E.g. 

2 



k A E\(k) 



a x (k)e ik - x -al(k)e- lk - x )dk , 



ik-x 



(78) 



recalling that h = c = 1. 

This notation, A A , B A , E A , with the superscript A, will be used in this appendix only. 

For the first two examples we define v Xd (k) = ((k A e\{k))- /(27r) 3 / 2 |/c|, so that i(Lj(x) 
L){x)) = Bf{x)/^ for j = 1,2,3. 



Example 1: Assume that < w(y) < 1 for all y. Then Lemma Rl implies 



H f >^- [ B A {x) 2 w(x)dx - / ™(y)dy 

87T J R 3 8TT Z J K 3 



(79) 



To verify the norm condition (|75| ) we note first that in view of ( |72[) it suffices to assume 
that w(y) = 1. Then, w(k) = (27r) 3 / 2 5(k). On the right side of (f7|) we may use the equality 
YljVx,j(,k)v\> t j(k) = (27t) _3 5a,a' (because ((A; A £\(k)/\k\)j are the three components of two or- 



thonormal vectors). Thus, ( [Tq) is not only satisfied, it is also an identity with this choice of w. 
Finally, ( |74[ ) is exactly ( [79] ) since /| fc i< A \ k\dk = 7rA 4 . 

Example 2: Take w(y) = C5{x — y), with x some fixed point in R 3 and where C > is some 
constant. Then w(k) = C(2tt)~ 3 / 2 exp(ik ■ x) and \w(k)\ = C(2tt)~ 3 / 2 . We take v\j(k) as in 
Example 1. The right side of condition ( |75| ) equals CX)j | Jiju<A /a(^)^a,j'(^) exp(ifc • x)d/c| 2 , 
and a simple variational argument then says that we should choose /a(&) = exp(— ik ■ x)v\(k) • V, 



where V is some fixed vector. By spherical symmetry we can take f\(k) = V\i(k), and ( |75| ) 
becomes 



1>C/ El^iWI 2 ^ = c|.-^3-^A 3 . 

i|t|<A j 3 (2vr) J 



(80) 



With C = 9tHA- 3 , and with 



EE / l^,iW| 2 ^ / w(y)dy = 2 • (2^)- 3 • (4tt/4)A 4 C , 

" " ./|ft|<A JK 3 

(0) becomes H f >\-C ■ (2tt)- 1 B a (x) 2 - (l/8)vr- 2 • A 4 • C, or 



ff, > ^L A - 3 B A (x) 2 - -A 



(81) 



This can also be used [jD], 14] with x being the electron coordinate (which is an operator, to be 
sure, but is one that commutes with the field operators). 
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Example 3: If we replace (k A E\(k))j /\k\ in v\ j(k) by (e\(k))j then everything goes through 



as before and we obtain © and @ with £ A (x) 2 in place of B A ( 



x 



,2 



Example 4: We now take v X j(k) = (e A (/c)) i (2vr)- 3 / 2 /|^|, so that {L j (x)+L*{x)) = Af(x)/V2Tr 
for j = 1,2,3. The analysis proceeds as in Example 2, except that the normalization condition ( p0| ) 
becomes 

1>C I ^ Ail (fc)p^ = c|--^-47rA 

7|fe|<A , 3 (27r) d 



which leads to C = 37r 2 A _1 . We also have 



A 3 

so that (fr^l ) becomes 



E E / 1*1 K#)l 2 ^ / ™(z/)dj/ = 2 • (2^r 3 • (4^/2)A 2 C , 

A * ^|fc|<A JR 3 



^ > ^A-i^ A (x) 2 - -A 



4 



(82) 



C Appendix: Spectral properties of the Dirac operators 



In this appendix we sketch a proof of the fact that the operators Di(A) commute in the sense 
that all their spectral projections commute. First we start with some remarks concerning the 
self-adjointness of 

Di(A) = eti ■ {-iVi + y/aAfaj) + m(3 . (83) 

The subscript after a. is a reminder that the matrix acts on the spinor associated with the i-th 
particle. It is not easy to characterize the domain for this operator, but it is certainly defined and 
symmetric on TC% := C^°(IR 3 ) ^finite , where ^finite denotes the vectors in Fock space with finitely 
many photons. This is a dense subset oiTC^. We shall show that Di(A) is essentially self-adjoint on 
ri° N . To prove this we resort to a version of Nelson's commutator theorem given in |22[, Theorem 
X.37. 

Define the operator 

JV 

u = l + Y J (-^i + m 2 ) + AH f + A 2 , (84) 

i=l 

Observe that ^2iLi(—\ + rn 2 ) acts as a multiplication operator on Fourier space and Hf acts on 
the n photon component <& n (xi, . . . , xn', t\, ■ ■ ■ , tn; ki, . . . , k n , Ai, . . . , A n ) by 

multiplication with X^r=i u; (^)- The domain ri° N is a domain of essential self-adjointness for v. 
Certainly, v > 1 as an operator. 
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We shall show that there exists a constant c such that for all ^ £ 7~L%, 

\\Di(A)V\\ < c\\W\\ . (85) 

Certainly, 

\\Di(A)y\\ < (*,(-A J + m 2 )^) 1/2 + (^,^(x i ) 2 ^) 1/2 (86) 
and by Example 4 in Appendix |b| 

A(x) 2 < ^AH f + —A 2 . (87) 

The estimate 

||A0A)*|| < c\\v$\\ (88) 

follows easily from this. 

Next, we show that there exists a constant d such that for all ^ 6 

| (A(^)*, f*) - (f*, A(>1)*) I < rf||^ 1/2 1'|| 2 • (89) 

Since — icti • Vj + mf3 commutes with v when applied to vectors in TL%, the above estimate reduces 
to 

| (on ■ A(xi)V, i/*) - (*/¥, a 8 • | < d\\v l / 2 mf , (90) 

where we have dropped the fine structure constant. Since v as well as A(x) preserve Ti% and are 
symmetric, we can rewrite the above inequality as 

\(*,[ ai -A(xi),v]*) \ <<% 1/2 *|| 2 . (91) 
The commutator is the sum of 

3 

[cti ■ A(x), -A] = on ■ {d j A{x)d j + d j d j A(x)) =: X , (92) 
i=i 

and the operator 

-ikcxi ■ E{xi) = [on ■ A(xi),AHf] , (93) 
where E{x) is the electric field (65). By Schwarz's inequality, 



| XV) | < (*, (V.4) 2 ^)) + -A*)) , (94) 

and 

A| (V,E(x)V)\ < A(||^|| 2 + {^,E(x) 2 ^)) . (95) 
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By Example 3 in Appendix ^, it follows that as quadratic forms 

E(x) 2 < ^A 3 H f + ±A 4 (96) 
(VA)(x) 2 < AA^ + IA 4 . (97) 

The last inequality does not appear in Appendix |B| exactly as stated, but it can be derived in 
precisely the same fashion as the one for the magnetic field displayed there. The estimates (p4|)- 
(0) yield (provided A > 1) 

| (Di(A)V, vV) - (*/¥, Di(A)V) | (98) 
< C (*, {A 3 H f + A 4 )^) + (¥, (-A + A)*) (99) 
< CA 2 (¥,*/¥) , (100) 

for some constant C, which is the desired estimate. 

Thus, the operator Di(A) is essentially self-adjoint on TL%. This operator, being a sum of 
two self-adjoint operators, U{ := — ia-i ■ Vj + mf3 and Vi := Qj • y/aA(xi), is naturally defined on 
V{Ui) n T>{Vi) and is symmetric there. Since, TL° N C V{Ui) n T>(Vi) we also know that Di(A) is 
essentially self-adjoint on T>(Ui) D T>(Vi). Thus, by Theorem VIII. 31 in |21] the Trotter product 
formula is valid, i.e. 

e itD t (A) =s _ Um T .( t / m \m ; ( 10 1) 

m— >oo 

where 

Ti(t/m) := jit/raWi jit/mW _ ^q 2 j 

Certainly, the operator e*'^ commutes with e lsUj and e ts ^' , and likewise e itVi commutes with e tsU: > 
and e lsV i for all j 7^ i, and hence Ti(t/m) m commutes with Tj(s/n) n for all s,t,m and n. 
We shall use this to show the following 

C.l. LEMMA. For any two real numbers s and t the unitary groups e itDi ^ and e isD ^ A ^ com- 
mute. Moreover, this implies that the spectral projections associated with Di(A) and Dj{A) com- 
mute. 

Proof. For * € H N 

^ e itD l (A) e isD :i (A)^ _ e i8Dj(A)jtDi(A)tyn Q03) 

= linw+oo \\Ti(t/m) m e isD ^ A ^ - e isD ^Ti{t/m) m ^\\ (104) 
= linwoo lim^oo \\T i {t/m) m T i {s/n) n ^ - T^s/nYT^/m)™^ = . (105) 

The statement about the spectral projections follows from Theorem VIII. 13 in [21]. □ 
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D Appendix: Projections and symmetries 

The difficulty in defining the physical space 7i^ JS comes from the fact that the projection onto 



the positive energy subspace acts also on the Fock space. This is in contrast to 17] where no such 
problem arises. There the action of the permutation group obviously commuted with the projection 
onto the positive energy subspace. In our more general setting the commutation is still true but 
an explanation is needed, which we try to give with a minimal amount of formality. 

First, consider the one particle space Jt\ = L 2 (IR 3 ;C 4 ) (^) J- . The Dirac operator, as shown in 
Appendix |C], is a self-adjoint operator on TL\ and we denote its projections onto the positive and 
negative energy subspace by P + and P~. Note that P + + P~ is the identity. As explained in 
Section [|, the two projections are unitarily equivalent via 

p- = U*P + U = -UP + U (106) 

where, as in 



U =[ J . (107) 




The projection onto the positive energy subspace associated with the Dirac operator Di(A) is 
defined in the following fashion. Consider the vector ^ as in ([l2] ) and fix N — 1 x's and t's, namely 
all those except X{ and Tj. For almost every such choice (with respect to Lebesgue measure) the 
vector ^ defines a vector in H\. We know how the P ± act on such a vector and the extension to 
TCn we denote by P^^- It was shown in Appendix |C| that these spectral projectors commute with 
each other. 

Other interesting operators on T~Ln are the permutations. A permutation Peri^, for example, 
just exchanges the electron labels 1 and 2. From what has been explained above we have the 
formula 

Per^Pf = jfPer^ . (108) 

An immediate consequence is that P ± := IlfL 1 P^ 1 commutes with permutations. 
From this it follows that ( |14| ) can be rewritten as 

H P x JS := AP + H N = P + AH N . (109) 

We now address the question whether Ti^ ys is trivial or not. Denote by K, the subspace of 
TCn which consists of antisymmetric vectors ^ with the property that Uj^S? = for j = 1, . . . , N . 
The operators Uj are defined in (|39[) . Certainly each Uj has eigenvalues i and —i. The space fC 
is certainly infinite dimensional. It contains, e.g., determinantal vectors in L 2 (M 3 ;C ) tensor the 
photon vacuum. 
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D.l. LEMMA (H phys * is large). The space H p ^ ys is infinite dimensional. 

Proof. We shall show that 2 N / 2 P+ is an isometry from K, into TCj^ ys . Let I be a subset of the 
integers {1, . . . , N} and let J be its complement. Let 

Pi = U ie iP-U jeJ P; . (110) 

Note that J2i Pi = identity. Note also that 

Pi = (-)^UiU ie iP t + UiU je jP j + (111) 

which implies that ||-Pr*|| = ||P + ^||. This shows in particular that 

HP+^II 2 = 2" Ar ||^|| 2 , (112) 

which proves the isometry. □ 

Since we always consider the symmetric operator H^ ys in the sense of quadratic forms, it is 
necessary to construct a domain, Qn that is dense in 1-L v ^ ys and on which every term in H^ ys 
has a finite expectation value. Once it is shown that the quadratic form associated with H^ ys is 
bounded below, it is closable and its closure defines a self-adjoint operator, the Friedrich's extension 
of H^ ys . 

We first start with a technical lemma that will allow us to approximate vectors in H p ^ ys . 
D.2. LEMMA. For any f with 

/oo 
|/(t)|(l + |t|)dt<oo (113) 
-oo 

and 

f(D(A)) = J e- itD Wf(t)dt (114) 

we have that 

\\y/l + Hff(D(A))V\\ < max{Vl + 9A/8, ^J^h?' 2 }C f \\^/l + H f V\\ , (115) 
for all * € Hi. 

Proof. We shall assume that ^ is normalized. Since 

\\y/T+H- f f(D(A))n < J \f(t)\\\V^+H-f^ UD{A) ndt , (116) 
it suffices to prove the estimate 

K{t) := \\y/T+H~ f e"'"^* || < C(l + \t\)\\^l + H~^\\ . (117) 
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A simple calculation yields 

^K\t) = (e- itD Wy, i[H f , D{A)]e- ttD ^^ = - (e- iW ^y, a ■ E(x)e~ itD ^^j . (118) 

Here E(x) is the electric field 

. 2 

E{x) = ^-J2 dke x (k)y^{k) (e ik -*a x (k) - e^a^k)) . (119) 

By Schwarz's inequality 



^-K 2 {t) < (e- iW My,E(x) 2 e- iW Mn>y /2 . (120) 



By Example 3 in Section [B] 



E( x f < —A 3 H f + -A 4 (121) 



and hence K (t) satisfies the differential inequality 

-^K 2 (t) < (AK 2 (t) + 5) 1/2 , (122) 
where A = ^A 3 and B = ^A 4 . This can be readily solved (using K(t) > 1) to yield the estimate 

K(t) < (1 + B/A) 1 / 2 K(0) + VAt . (123) 

Thus 

\\ y /T+H} e - itD ( A >V\\ < C(l + liDHv/r+fl/^H , (124) 
where C is the maximum of (1 + B/A) l / 2 and \[~A. 

□ 

Next we consider a sequence of functions f n G C^°((0,oo)) everywhere less or equal to 1, such 
that f n is identically equals to 1 on the interval [l/n,n]. Clearly, as n — > oo, Tif =l f n (Di{A)) —>■ 
P + strongly in TLn and hence AH^ =1 f n (Di(A)) — » AP + strongly in TLn- We denote the range 
of AT\.f =l fn{Di{A)) restricted to the subspace of TLn consisting of states with finite field energy 
expectation by Q N . Finally we define the domain Qn = U^L 1 Q^ r . Together with Lemma D.2 we 
have the following Corollary. 

D.3. COROLLARY. The domain Qn is dense in ?^ ys . Moreover for any vector ^> G Qn the 
field energy Hf has finite expectation value. 

Proof. Simply note that the functions /„ have a rapidly decaying Fourier transform for each n. 
Therefore, by Lemma D.2 the field energy has a finite expectation value for any vector ^ £ Q N . 
Note, as before, the antisymmetrization operator A commutes with U^L 1 f n (Di(A)). Thus, the 
field energy has finite expectation value for any ^ 6 Qn- The density of Qn in TL.ffi*' was shown 
before. □ 
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Now we are ready to state the main lemma of this section. 

D.4. LEMMA. For every £ Qn, the Dime operators Di(A), the Coulomb potential V c and the 
field energy Hf have finite expectation values. Thus, H^ ys is defined as a quadratic form on Qn 
which is dense in W^ iys . 

Proof. The operators Di(A) 2 have finite expectation values on Qn- They are of the form T P (A) = 
[(p + yfaA(x)) 2 + yjao ■ B(x)] (g) h where I2 is the 2x2 identity. By (|8l|) the magnetic field is 
bounded by the field energy and hence has finite expectation value on the domain Qn- Thus, 
+ TfaA(xi)) 2 (g> I2 also has finite expectation value on Qn for i = 1,... ,N, and hence the 
Coulomb potential V c , which is relatively bounded with respect to X^Li(p» + ^faA(xi)) 2 has finite 
expectation values on Qn- □ 



E Appendix: Various forms of instability 

In the introduction we talked about the need of using the positive spectral subspace of the Dirac 
operator D(A), which includes the magnetic vector potential; this led to all sorts of complications 



in the analysis leading to our main stability Theorem 2.1. In this section we show that various 
models in which an electron is defined, instead, by the positive spectral subspace of the free Dirac 
operator -D(O) are unstable. In the case of a classical magnetic field such an analysis was carried 
out in |l7]] and greatly simplified in ||. Also, in Q the stability analysis was carried out for a 
quantized radiation field without a cutoff. In what follows, we rely mostly on the work in ||. We 
also show that the D(A) choice is unstable if Za or a is too large — as expected. 

All the results about stability and instability are summarized in the two tables in Section [l]. 
We remind the reader that instability of the first kind means that the Hamiltonian is unbounded 
below, while instability of the second kind means that it is bounded below but not by a constant 
times N + K. 



E.l Instability without Coulomb potential (free Dirac operator) 

Already the free problem, i.e., without Coulomb interactions, shows signs of instability. The Hamil- 
tonian is given by 

N 

H N = Y, D M) + H f ■ (125) 
j'=i 

If the field is classical Hf has to be replaced by (1/8tv) L 3 |S(x)| 2 dx as in (||). 

We consider first the case where the magnetic vector potential is classical. In particular the 
Hilbert space W^ ee is the antisymmetric tensor product of N copies of P + L 2 (1R 3 ; C 4 ), i.e., the part 
of L 2 (M 3 ;C 4 ) that is in the positive spectral subspace of the free Dirac operator. Note that there 
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is no Fock space in this case. In ||] Theorems 1 and 3 the authors construct, for any N, a trial 
Slater determinant ip in Tl^ e , and a classical field A so that the energy is bounded above by 

(V, ## assica V) =: Sty, A) < aiV 4 / 3 - abN 2 , (126) 

where a and b are constants independent of N. The scaling 

ip -»■ ^ , and 4 -» A" (127) 

where 

• • • ,x N ;n, ...,t n ) = [i 3N ' 2 ip(nxx, . . .,fiX N ;Ti, ...,t n ) (128) 

and 

A»(x) = fiAiptx) (129) 

can be used to get the upper bound 

Sty» A 1 ") < pL (AN 1 '* - aBN 2 ^j . (130) 

Thus, by choosing N > aa~ 3 ^ 2 /b, the ground state energy is negative and can be driven to — oo by 
letting [j, — > oo, i.e., stability of the first kind is violated. 

Using coherent states, it was shown in || that this same result extends to the problem with a 
quantized magnetic vector potential without ultraviolet cutoff. 

If A(x) carries an ultraviolet cutoff the fi scaling argument cannot be applied. The energy, 
however, is not bounded below by const, x N, as we see from ( |126|) , and hence stability of the 
second kind is also violated — and this no matter how small a may be and whether the magnetic 
vector potential is quantized or not. The reader might wonder how to construct an A(x) that 
satisfies the conditions in || and, at the same time has an ultraviolet cut off A. Remark 1 on page 
1782 of H explains that almost any cutoff A(x) will suffice for the purpose. 

Notice, that when we use the positive spectral subspace of D{A) instead, the stability of the 
problem without Coulomb potential is completely trivial, since the Hamiltonian is positive, by 
definition. 



E.2 Instability with Coulomb potential (free Dirac operator) 

Adding the Coulomb potential complicates the analysis owing to the repulsion between the electrons 
which is present even if there are no nuclei. To some extent this positive energy is balanced by the 
electron- nuclei attraction if sufficiently many nuclei with sufficiently strong charges are present. It 
is shown in [17] that if 

K K 

Zj > (const.)cT 3 / 2 and ^ Z 2 > 2 , (131) 
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then the positions of the nuclei can be chosen such that the total Coulomb energy is negative. 
Thus, if in addition, iVa 3 / 2 is sufficiently large, stability of the first kind does not hold for classical 
magnetic vector potentials as well as for a quantized magnetic vector potential (without ultraviolet 
cutoff) — no matter how small a may be. 

The situation is more complicated when the field carries an ultraviolet cutoff. The main reason 
is that the field variable is no longer an active participant for driving the energy towards minus 
infinity, but it is an active participant in destroying stability of the second kind. We have 

E.l. LEMMA. Let a > and assume that \131 ) holds. Then the system using the projection 
onto the positive subspace of the free Dirac operator -D(O) is unstable of the second kind, even with 
an ultraviolet cutoff. 

Proof. The lemma follows immediately from (13C) and ( |131| ) together with the observation in _] 
on how to use coherent states to carry these results over to the quantized field case. □ 



Lemma E.l is the main reason that the restriction to the positive spectral subspace of D(0) is 
inadequate for a model of matter interacting with radiation. 

The main result of this paper is the stability of the second kind for the system where the positive 
subspace of the Dirac operator D{A) is used. This result holds provided that maxj ZjO. and a is 
sufficiently small. Our final Lemma |E.2| shows that that the two conditions max,,- Zja small and a 
small are in fact necessary. It suffices to show this for the case of one electron interacting with K 
nuclei, each having charge Z. We have to assume that the electron mass is strictly positive, but we 
suspect that this assumption is technical and not needed. 

E.2. LEMMA. Assume that Za > 4/7r and m > 0. Then the one-electron Hamiltonian 



^Phys = p+ . D{A) _ Za ^ _____ + z , a ^ + Hf X P+ (132) 

\ 3=1 3 i<j J / 

is unbounded below. Moreover, there is a number a c such that for any fixed a > a c and any fixed 
Z > 0, there exists K sufficiently large so that this Hamiltonian is unbounded below. 

For the classical instead of the quantized A field see || ___\, _], |24|, ___] . 

Proof. The idea is to reduce this problem to the relativistic one without spin and without radiation 
field. Set D + = P+D(A) and L>_ = -P~D(A), so that D(A) = D + -D_ and \D(A)\ =D + + D_. 
As a trial function we pick tp = g<3 10), where g is a spinor that satisfies Ug = ig, recalling the defini- 
tion of U from Appendix [D| ( |I06| , |I07[ ), and where |0) is the photon vacuum. A straightforward cal- 
culation (which repeatedly uses Schwarz's inequality and the facts that P~ = U*P + U = —UP + U, 
and hence P~ip = —i UP + ip) shows that 

K 



D + - Za Y_] I p| + Z 2 a ^ — — + H f 

j = l J Ky J 1 



p+V> 
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2 I P V>, 



+ 2' P ^' 



+ (^,P + H f P + ^) 



< 2l P ^ 



if 
K 



1 




- R 1 




1 




- i? 1 




1 


\x 


- 




1 


\x 


- Rj\ 



i<j 



\Ri 


- RA 




l 


\Ri 


- RA 




1 


\Ri 


- Rj\ 




1 



p + 1p 



P'tp 



j=l J i<j 

2 ^-f b-i?,- ^ Ri-RA 

3=1 J i<j ' 

Re f P+V>, ^ £ J^J^P'^j + (V>, P + H f P + i>) 



P+ip 



< ^ IV, 



v ;| 9 Z^ u _ p J Z^ p, _ 



2 ^ s - i? 7 l 



^ + (^,P + H f P + 4>) 



The lemma will be proved by showing that the first term in the last expression can be made 
as negative as we like while the field energy term is uniformly bounded. Note, as in (^) that 
|D(A)| = \JT P (A) + m 2 . The operator inequality 



(p + V^A(x)) 2 < (1 + e)p 2 + (1 + ^)aA(xf 



follows easily from Schwarz's inequality, for any e > 0. From (|8^) we have that 

2 



and from (pi) we have that 



A(xr < —AH f + -A 2 
Sir 7r 



B {xf < —A 3 H f + -A 4 . 

V 1 ~ 97T ; 7T 



Using the operator monotonicity of the square root it follows that 



T P (A) + m 2 < J(l + e)p 2 + X £ AH f + Y £ A 2 + m 2 

Where X £ and Y £ are constants that tend to infinity as e tends to zero. 
Thus, recalling that ip = g ® \0) 

(V, \D(A)\i/i) < (g, + e)p 2 + Y £ A 2 g 

The remaining task is to analyze the quadratic form 



(133) 



(134) 



(135) 



(136) 



V(l + e)p 2 + Y £ A 2 + m 2 



Za 



K 1 1 
Z^ u _ p, Z^ p, _ 



(137) 



(138) 
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For any fixed e, A and m the terms Y £ A 2 + m 2 can be scaled away, and this leads to the quadratic 
form 



\Ri — Rj 



(139) 

j=l J Kj 

which has been analyzed in detail. Kato |Tl]] showed that instability of the first kind occurs if 
Za/2y/\ + e > 2/it which yields our first stated condition for instability. Later on, it was shown in 
H that there exists a c so that for a > a c and for any Z > there exists K so that an instability 



of the first kind occurs. See also [19], Theorem 3. 

Next we address the field energy term and it is here where the assumption about the positive 
mass comes in. The projection P + can be written as P + = (1/2) (/ + D(A)/\D(A)\) and, since ip 
is proportional to the vacuum, 



(ip,P + H f P + ip) 



1 (.,. D ( A ) „ D ( A ) 



H 



D(A)r f '\D(A)\ 



a\(k), 



\D(A)\ 



dk . 



With the help of the expression 



D(A) 



D(A) 1 



\D(A)\ 7rJ t + D(Af^t 



dt , 



(140) 



and the fact that (t + D(A) 



2^-l _ l 



(D(A) - iVi)" 1 + \ {D(A) + iVly 1 , the commutator in the 



last expression of (|E.2| ) can be written as 



a x (k). 



\D(A)\ 



2vr J D(A) - iy/i 



2vr J d(A) + iy/i 



[a x (k),D(A)] 
[a x (k),D(A)] 



1 



1 



D(A) -ixftyft 

1 1 

D{A) + iViVi 



dt 
dt . 



(141) 



and 



[a x (k),D(A)] 



e x (k) t 

y/^ik)' 



-ik-x 



where xa(^) denotes the characteristic function of the ball in k space that has radius A and is 
centered at the origin. Hence, 



D(A) 



\D(A)\ 
and the estimate 

follows. 



< 



1 1 



1 1 



y/uAk) vr J m 2 +ty/t 



(i/;,P + H f P + ip) < 



dtU\\xA(k) 



1 



m 



4vrA 3 
3m 2 



XA(k) , (142) 

(143) 
□ 
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